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BACKGROUND OF THE INVENTION 
1. FIELD OF THE INVENTION 

The present invention is generally in the field of speech coding. In particular, the 
present invention is in the field of noise suppression for speech coding purposes. 
5 2. BACKGROUND ART 

Today, noise reduction has become the subject of many research projects in 
various technical fields. In the recent years, due the tremendous demand and growth in 
2 the areas of digital telephony, the Internet and cellular telephones, there has been an 
Iji intense focus on the quality of audio signals, especially reduction of noise in speech 
1# signals. The goal of an ideal noise suppressor system or method is to reduce the noise 
T level without distorting the speech signal, and in effect, reduce the stress on the listener 
CO and increase intelligibility of the speech signal. 

f« Technically, there are many different ways to perform the noise reduction. One 

noise reduction technique that has gained ground among the experts in the field is a noise 

15 reduction system based on the principles of spectral weighting. Spectral weighting means 
that different spectral regions of the mixed signal of speech and noise are attenuated or 
modified with different gain factors. The goal is to achieve a speech signal that contains 
less noise than the original speech signal. At the same time, however, the speech quality 
must remain substantially intact with a minimal distortion of the original speech. Another 

20 important design consideration is that the residual noise, i.e. the noise remaining in the 
processed signal, must not sound unnatural. 

Typically, the spectral weighting technique is performed in the frequency domain 
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using the well-known Fourier transform. To explain the principles of spectral weighting 
in simple terms, a clean speech signal is denoted with s(k), a noise signal is denoted with 
n(k), and an original speech signal is denoted with o(k), which may be formulated as o(k) 
= s(k) + n(k). Now, taking the Fourier transform of this equation leads to 0(f) = S(f) + 

5 N(£). At this step, the actual spectral weighting may be performed by multiplying the 
spectrum 0(f) with a real weighting function, such as W(f) >= 0. As a result, P(f) = W(f) 
0(f), and the processed signal p(k) is obtained by transforming P(f) back into the time 

•4 domain. Now, below, a more elaborate system 100, including a conventional noise 

P.' suppression module 106 is discussed. The conventional noise suppression module 106 of 

m the speech pre-processing system 100 is that of the Telecommunication Industry 
Association Interim Standard 127 ("IS-127"), which is known as Enhanced Variable Rate 

i j Coder ("EVRC"). The IS-127 specification is hereby fully incorporated by reference in 

[P the present application. 

As stated above, FIG. la illustrates a conventional speech pre-processing system 

15 100, which includes a noise suppression module 106. After reading and buffering 
samples of the input speech 101 for a given speech frame, an input speech signal 101 
enters the speech preprocessor system 100. The input speech signal 101 samples are then 
analyzed by a silence enhancement module 102 to determine whether the speech frame is 
pure silence, in other words, whether only silence noise is present. Next, the silence 

20 enhanced input speech signal 103 is scaled down by the high-pass filter module 104 to 
condition the input speech 101 against excessive lose frequency that degrade the voice 
quality. 
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The high-pass filtered speech signal 105 is then routed to a noise suppression 
module 106. The noise suppression module 106 performs a noise attenuation of the 
environmental noise in order to improve the estimation of speech parameters. 

The noise suppression module 106 performs noise processing in frequency domain 
5 by adjusting the level of the frequency response of each frequency band that results in 
substantial reduction in background noise. The noise suppression module 106 is aimed at 
improving the signal-to-noise ratio ("SNR") of the input speech signal 101 prior to the 
% speech encoding process. Although the speech frame size is 20 ms, the noise suppression 
l;| module 106 frame size is 10 ms. Therefore, the following procedures must be executed 
tq: two times per 20 ms speech frame. For the purpose of the following description, the 
J current 10 ms frame of the high-pass filtered speech signal 105 is denoted m. 
W As shown, the high-pass filtered speech signal 105, denoted {S hp (n)}, enters the 

n first stage of the noise suppression module 106, i.e. Frequency Domain Conversion stage 

if"! 

110. At the frequency domain conversion stage 110, S hp (n) is windowed using a 
15 smoothed trapezoid window, in which the first D samples of the input frame buffer 
{d(m)} are overlapped from the last D samples of the previous frame, where this overlap 
is described as: d(m,n) = d(m-l,L+n); 0 < n < D, where m is the current frame, n is the 
sample index to the buffer {d(m)}, L = 80 is the frame length, and D = 24 is the overlap 
or delay in samples. The remaining samples of the input buffer {d(m)} are then pre- 
20 emphasized at the Frequency Domain Conversion stage 110 to increase the high to low 
frequency ratio with a pre-emphasized factor Cp = -0.8 according to the following: 
d(m,D+n) = S kp (n) + C P S hp (n-1); 0 < n < L. This results in the input buffer containing L 
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+ D = 104 samples in which the first D samples are the pre-emphasized overlap from the 
previous frame, and the following L samples are pre-emphasized input from the current 
frame m. 

Next, a smoothed trapezoidal window is applied to the input buffer (d(m)} to form 
5 a Discrete Fourier Transform ("DFT") data buffer (g(n)}, defined as: 



d(m,n) sin (ic(n+0.5)/2D) ;0<n<D, 

d(m,n) ;D<n<L, 

d(m,n) sin 2 (7r(n-L+D+0. 5)/2D) ;0<n<D, 

0 ; D+L <n<M, 



f Z where M =128 is the DFT sequence length. At this point, a transformation of g(n) to the 
p frequency domain is performed using the DFT to obtain G(k). A transformation 
5 technique, such as a 64-point complex Fast Fourier Transform ("FFT") may be used to 
|| convert the time domain data buffer g(n) to the frequency domain data buffer spectrum 
O G(k). Thereafter, G(k) is used to computer noise reduction parameters for the remaining 
blocks, as explained below. 

The frequency domain data buffer spectrum G(k) resulting from the Frequency 
Domain Conversion stage 110 is used to estimate channel energy E ch (m) for the current 
20 frame m at Channel Energy Estimator stage 115. At this stage, the 64-point energy bands 
are computed from the FFT results of stage 101, and are quantized into 16 bands (or 
channels). The quantization is used to combine low, mid, and high frequency 
components and to simplify the internal computation of the algorithm. Also, in order to 
maintain accuracy, the quantization uses a small step size for low frequency ranges, 
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increased the step size for higher frequencies, and uses the highest step size for the 
highest frequency ranges. 

Next, at Channel SNR Estimator stage 120, quantized 16-channel SNR indices 
a q (i) are estimated using the channel energy E ch (m) from the Channel Energy Estimator 
5 stage 115, and current channel noise energy estimate E n (m) from Background Noise 
Estimator 140 which continuously tracks the input spectrum G(k). In order to avoid 
undervaluing and overvaluing of the SNR, the final SNR result is also quantized at the 
hJ Channel SNR Estimator 120. Then, a sum of voice metrics v(m) at Voice Metric 
m Calculation stage 130 is determined based upon the estimated quantized channel SNR 
W indices a q (i) from the Channel SNR Estimator stage 120. This process involves a 
J transformation of the actual sum of all sixteen signal-to-noise ratios from a predetermined 
M voice metric table with the quantized channel SNR indices a q (i). The higher the SNR, the 
f{ higher the voice metric sum v(m). Because the value of the voice metric v(m) is also 

quantized, the maximum and the minimum values are always ascertainable. 
15 Thereafter, at Spectral Deviation Estimator stage 125, changes from speech to 

noise and vice versa are detected which can be used to indicate the presence of speech 
activity of a noise frame. In particular, a log power spectrum E db (mJ) is estimated based 
upon the estimated channel energy E ch (m), from the Channel Energy Estimator stage 115, 
for each of the sixteen channels. Then, an estimated spectral deviation A E (m) between a 
20 current frame power spectrum E db (m) and an average long-term power spectral estimate 
E db (m) is determined. The estimated spectral deviation A E (m) is simply a sum of the 
difference between the current frame power spectrum E db (m) and the average long-term 
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power spectral estimate E db (m) at each of the sixteen channels. In addition, a total 
channel energy estimate E tot (m) for the current frame is determined by taking the 
logarithm of the sum of the estimated channel energy E ch (m) at each frame. Thereafter, 
an exponential windowing factor a(m) as a function of the total channel energy E tot (m) is 

5 determined, and the result of that determination is limited to a range determined by a 
predetermined upper and lower limits a H and a L9 respectively. Then, an average long- 
term power spectral estimate for the subsequent frame E db (m+l f i) is updated using the 

U exponential windowing factor a(m) 9 the log power spectrum E db (m) 9 and the average long- 

y=j term power spectral estimate for the current frame E db (m). 

W~ With the above variables determined at the Spectral Deviation Estimator stage 125, 

* ?s noise estimate is updated at Noise Update Decision stage 135. At this stage 135, a noise 
Co frame indicator update Jlag indicating the presence of a noise frame can be determined 
::{ by utilizing the voice metrics v(m) from the Voice Metric Calculation stage 130, and the 
total channel energy E tot (m) and the spectral deviation A E (m) from the Spectral Deviation 
15 Estimator stage 125. Using these three pre-computed values coupled with a simple delay 
decision mechanism, the noise frame indicator update Jlag is ascertained. The delay 
decision is implemented using counters and a hysterisis process to avoid any sudden 
changes in the noise to non-noise frame detection. The pseudo-code demonstrating the 
logic for updating the noise estimate is set forth in the above-incorporated IS- 127 
20 specification and shown in FIG. lb. 

Now, having updated the background noise at the Noise Update Decision stage 
135, at Channel Gain Calculation stage 150, it is determined whether channel SNR 
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modification is necessary and whether to modify the appropriate channel SNR indices 
a q (i). In some instances, it is necessary to modify the SNR value to avoid classifying a 
noise frame as speech. This error may stem from distorted frequency spectrum. By 
analyzing the mid and high frequency bands at Channel SNR Modifier stage 145, the pre- 
5 computed SNR can be modified if it is determined that a high probability of error exists in 
the processed signal. This process is set forth in the above-incorporated IS- 127 
specification, as shown in FIG. lc. 
q Referring to FIG. lc, the quantized channel SNR indices a q (i) determined at the 

Y 1 Channel SNR Estimator 120 are verified to be greater or equal to a predetermined channel 
lb] SNR index threshold level, i.e. INDEXJTHLD, which is set at 12. Thereafter, if it is 
I determined that the index counter is less than a predetermined index counter threshold 
Jjt level (INDEX_CNT_THLD=5), a channel SNR modification flag may be set to indicate 
Q that the channel SNR must be modified and the channel SNR indices o q (i) are modified to 
obtain modified channel SNR indices a' q (i) or the channel SNR modification flag may be 
15 reset to indicate that the modification is not necessary, and the modified channel SNR 
indices a q (i) are not changed from the original values a f q (i) = a q (i). 

Now, if the voice metric sum v(m) determined at the Voice Metric Calculation 
stage 130 is determined to be less than or equal to a predetermined metric threshold level, 
i.e. METRIC _THLD=4 5, or if the channel SNR indices a q (i) are less than or equal to a 
20 predetermined setback threshold level, i.e. SETS A CKJTHLD =12, the modified channel 
SNR indices a f q (i) are set to one. Else, the modified channel SNR indices a f q (i) are not 
changed from the original values a f q (i) = a q (i). In the following segment, in order to limit 
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the modified channel SNR indices a q (i) to an SNR threshold level a th , it is first 
determined whether the modified channel SNR indices o' q (i) are less than the SNR 
threshold level o th . If so, the threshold limited and modified channel SNR o" q (i) indices 
are set to the threshold level a th , i.e. (a" q (i) = o i} ). Else, the SNR indices a" q (i) are not 
5 changed, i.e., <r" q (i) = a\(i). 

Turning back to FIG. la, the threshold limited, modified channel SNR indices 
a" q (i) are provided to the Channel Gain Calculation stage 150 to determine an overall gain 
*\i factor y n for the current frame based upon a pre-set minimum overall gain y mm a noise 
Y 1 floor energy E floor , and the estimated noise spectrum of the previous frame E n (m-1). Next, 
M}j the channel gain in the db domain, i.e. y^(i), is determined based on the following 
* Si equation: 



□ where the gain slope ju g is constant factor, set to 0.39. In the following stage, the channel 
gain * s converted from the db domain to linear channel gains, i.e. y c *(i), by taking 



15 the inverse logarithm of base 10, i.e. y ch (i) = minj^ 1, io ydb(!)/2 ° j>. Therefore, for a given 
channel, y ch has a value less than or equal to one, but greater than zero, i.e. 0 < y C hO) - 1- 
The gain y ch should be higher or closer to 1 .0 to preserve the speech quality for strong 
voiced areas and, on the other hand, the gain y ch should be lower or closer to zero to 
suppress noise in noisy areas. Next, the linear channel gains y ch (i) are applied to the G(k) 

20 signal by a gain modifier 155 producing a noise-reduced signal spectrum H(k). Finally, 
H(k) signal is converted back into time domain at Time Domain Conversion stage 160 
resulting in noise reduced signal S f (n) in the time domain. 



y db (i) = p g ( a" q (i) - a th ) +y n ;0<i<N c 
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The above-described conventional approach, however, is a simplistic approach to 
noise suppression, which only considers one dynamic parameter, i.e. the dynamic change 
in the SNR value, in determining the channel gains y ch (i). This simplistic approach 
introduces various downfalls, which may in turn cause a degradation in the perceptual 
5 quality of the voice signal that is more audible than the noise signal. The shortcomings 
and inaccuracies of the conventional system 100, which are due to its sole reliance on the 
SNR value, stem from the facts that the SNR calculation is merely an estimation of the 
U noise to signal, and that the SNR value is only an average, which by definition may be 
)~:\ more or less than the true SNR value for specific areas of each channel. As a result of its 
Vf mere reliance on the SNR value, the conventional approach suffers from improperly 
J altering the voiced areas of the speech, and thus, causes degradation in the voice quality. 
K Accordingly, there is an intense need in the art for a new and improved approach 

)ll to noise suppression that can overcome the shortcomings in the conventional approach 
and produce a noise-reduced speech signal with a superior voice quality. 

15 
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SUMMARY OF THE INVENTION 
In accordance with the purpose of the present invention as broadly described 
herein, there is provided method and system for suppressing noise to enhance signal 
quality. 

5 According to one aspect of the present invention, an input signal enters a noise 

suppression system in a time domain and is converted to a frequency domain. The noise 
suppression system then estimates a signal to noise ratio of the frequency domain signal. 
J;f Next, a signal gain is calculated based on the estimated signal to noise ratio and a voicing 
|j1 parameter. In one aspect of the present invention, the voicing parameter may be 
determined based on the frequency domain signal. 

In another aspect, the voicing parameter may be determined based on a signal 
W ahead of the frequency domain signal with respect to time. In that event, the voicing 
S parameter is fed back to the noise suppression system to calculate the signal gain. 

After calculating the gain, the noise suppression system modifies the signal using 
15 the gain to enhance the signal quality. In one aspect, the modified signal may be 
converted from the frequency domain to time domain for speech coding. 

In one aspect, the voicing parameter may be a speech classification. In another 
aspect, the voicing parameter may be a signal pitch information. Yet, the voicing 
parameter may be a combination of several speech parameters or a plurality of parameters 
20 may be used for calculating the gain. In yet another aspect, the voicing parameter(s) may 
be determined by a speech coder. 

In one aspect of the present invention, the signal gain may be calculated based on 
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7db = Mg ( °"q ~ °th) + 7m such that ju g is adjusted according to the voicing parameter(s). In 
other aspects, the voicing parameter(s) may be used to adjust other parameters in the 
above-shown equation, such as a th or y„,or elements of any other equation used for noise 
suppresion purposes. 

Other aspects of the present invention will become apparent with further reference 
to the drawings and specification, which follow. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The features and advantages of the present invention will become more readily 
apparent to those ordinarily skilled in the art after reviewing the following detailed 
description and accompanying drawings, wherein: 

FIG. la illustrates a conventional speech pre-processing system; 

FIG. lb illustrates a conventional pseudo-code for implementing the Noise Update 
Decision module of FIG. la; 

FIG. lc illustrates a conventional pseudo-code for implementing the Channel SNR 
Modifier module of FIG. la; 

FIG. 2 illustrates a speech processing system according to one embodiment of the 
present invention; 

FIG. 3 illustrates voiced, unvoiced and onset areas of a speech signal in time 
domain; and 

FIG. 4 illustrates a speech signal in frequency domain. 
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DETAILED DESCRIPTION OF THE INVENTION 
The present invention discloses an improved noise suppression system and 
method. The following description contains specific information pertaining to the 
Extended Code Excited Linear Prediction Technique ("eX-CELP"). However, one 
skilled in the art will recognize that the present invention may be practiced in conjunction 
with various speech coding algorithms different from those specifically discussed in the 
present application. Moreover, some of the specific details, which are within the 
knowledge of a person of ordinary skill in the art, are not discussed to avoid obscuring the 
present invention. 

The drawings in the present application and their accompanying detailed 
description are directed to merely example embodiments of the invention. To maintain 
brevity, other embodiments of the invention which use the principles of the present 
invention are not specifically described in the present application and are not specifically 
illustrated by the present drawings. 

FIG. 2 illustrates a block diagram of an example encoder 200 capable of 
embodying the present invention. As shown, the encoder 200 is divided into a speech 
pre-processor block 210 and a speech processor block 250. First, an input speech signal 
201 enters the speech pre-processor block 210. After reading and buffering samples of 
the input speech 201 for a given speech frame, the input speech signal 201 samples are 
analyzed by a silence enhancement module 202 to determine whether the speech frame is 
pure silence, in other words, whether only silence noise is present. 

The silence enhancement module 202 adaptively tracks the minimum resolution 
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and levels of the signal around zero. According to such tracking information, the silence 
enhancement module 202 adaptively detects, on a frame-by-frame basis, whether the 
current frame is silence and whether the component is purely silence noise. If the silence 
enhancement module 202 detects silence noise, the silence enhancement module 202 
5 ramps the input speech signal 201 to the zero-level of the input speech signal 201. 
Otherwise, the input speech signal 201 is not modified. It should be noted that the zero- 
level of the input speech signal 201 may depend on the processing prior to reaching the 
y encoder 200. In general, the silence enhancement module 202 modifies the signal if the 
iv\ sample values for a given frame are within two quantization levels of the zero-level. 

In short, the silence enhancement module 202 cleans up the silence parts of the 

k . if 

; " input speech signal 201 for very low noise levels and, therefore, enhances the perceptual 
ra quality of the input speech signal 201 . The effect of the silence enhancement module 202 
U becomes especially noticeable when the input signal 201 originates from an A-law source 
or, in other words, the input signal 201 has passed through A-law encoding and decoding 
1 5 immediately prior to reaching the encoder 200. 

Continuing with FIG. 2, the silence enhanced input speech signal 203 is then 
passed through a high-pass filter module 204 of a 2 nd order pole-zero with a cut-off 
frequency of 240 Hz. The silence enhanced input speech signal 203 is scaled down by a 
factor of two by the high-pass filter module 204 that is defined by the following transfer 
20 function. 

0.92727435-1.8544941z~' +0.92727435z' 2 
1 - 1 .9059465z + 0.9 1 14024z ~ 2 
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The high-pass filtered speech signal 205 is then routed to a noise suppression 
module 206. At this point, the noise suppression module 206 attenuates the speech signal 
in order to provide the listener with a clear sensation of the environment. As shown in 
FIG. 2, the noise suppression module 206, including a channel gain calculation module 
5 208 receives a number of voicing parameters from the speech processor block 250 via a 
voicing parameter feedback path 260. The voicing parameters include various speech 
signal parameters, such as speech classification, pitch information, or any other 
L \ parameters that are calculated by the speech processor block 250 while processing the 
Lri input speech signal 201. The voicing parameters are then fed back into channel gain 
K? calculation module 208 of the noise suppression module 201 to compute the gain {y c hQ)}-> 
I so as to improve the speech quality. This process is discussed in more details below. 
f fi Next, as the pre-processed speech signal 207 emerges from the speech pre- 

S processor block 210, the speech processor block 250 starts the coding process of the pre- 
processed speech signal 207 at 20ms intervals. At this stage, for each speech frame 
15 several parameters are extracted from the pre-processed speech signal 207. Some 
parameters, such as spectrum and initial pitch estimate parameters may later be used in 
the coding scheme. However, other parameters, such as maximal sample in a frame, zero 
crossing rates, LPC gain or signal sharpness parameters may only be used for 
classification and rate determination purposes. 
20 As further shown in FIG. 2, the pre-processed speech signal 207 enters a linear 

predictive coding ("LPC") analysis module 220. A linear predictor is used to estimate the 
value of the next sample of a signal, based upon a linear combination of the most recent 
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sample values. At the LPC analysis module 220, a 10 th order LPC analysis is performed 
three times for each frame using three different-shape windows. The LPC analyses are 
centered and performed at the middle third, the last third and the look-ahead of each 
speech frame. The LPC analysis for the look-ahead is recycled for the next frame as the 
LPC analysis is centered at the first third of each frame. Accordingly, for each speech 
frame, four sets of LPC parameters are available. 

A symmetric Hamming window is used for the LPC analyses of the middle and 
last third of the frame, and an asymmetric Hamming window is used for the LPC analysis 
of the look-ahead in order to center the weight appropriately. For each of the windowed 

AM 

segments the 10 order, auto-correlation is calculated according to r(k) ^^s^n) • s w (n- k) , 



where s w (n) is the speech signal after weighting with the proper Hamming window. 

Bandwidth expansion of 60Hz and a white noise correction factor of 1.0001, i.e. 
adding a noise floor of -40dB, are applied by weighting the auto-correlation coefficients 
according to r w (k) = w(k) • r(k) , where the weigthing function is given by 



Based on the weighted auto-correlation coefficients, the short-term LP filter 



coefficients, i.e. A(z) = \-^a J *z ~' , are estimated using the Leroux-Gueguen algorithm, 



and the line spectrum frequency ("LSF") parameters are derived from the polynomial 
A(z). The three sets of LSFs are denoted lsf y <l), it = i,2...,io, where lsf 2 (£), lsf 3 (A), and 



1.0001 



* = 0 




10 
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lsf 4 (£)are the LSFs for the middle third, last third and lookahead of each frame, 
respectively. 

Next, at the LSF smoothing module 222, the LSFs are smoothed to reduce 
unwanted fluctuations in the spectral envelope of the LPC synthesis filter (not shown) in 

5 the LPC analysis module 220. The smoothing process is controlled by the information 
received from the voice activity detection ("VAD") module 224 and the evolution of the 
spectral envelope. The VAD module 224 performs the voice activity detection algorithm 

^ for the encoder 200 in order to gather information on the characteristics of the input 

ti speech signal 201. In fact, the information gathered by the VAD module 224 is used to 

\4 control several functions of the encoder 200, such as estimation of signal to noise ratio 
("SNR"), pitch estimation, classification, spectral smoothing, energy smoothing and gain 

hi normalization. Further, the voice activity detection algorithm of the VAD module 224 
may be based on parameters such as the absolute maximum of frame, reflection 
coefficients, prediction error, LSF vector, the 10 th order auto-correlation, recent pitch lags 

1 5 and recent pitch gains. 

Continuing with FIG. 2, an LSF quantization module 226 is responsible for 
quantizing the 10 th order LPC model given by the smoothed LSFs, described above, in the 
LSF domain. A three-stage switched MA predictive vector quantization scheme may be 
used to quantize the ten (10) dimensional LSF vector. The input LSF vector (unquantized 

20 vector) originates from the LPC analysis centered at the last third of the frame. The error 
criterion of the quantization is a WMSE (Weighted Mean Squared Error) measure, where 
the weighting is a function of the LPC magnitude spectrum. The objective of the 
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quantization is set forth as {lsf n (i),isf n (i),...,isf n (iO)}=ar g min|£w 1 -(isf n (*)-isf n (*)) 2 j, where the 

weighting is w, =|/>(lsf„(0)|°\ where \p(f)\ is the LPC power spectrum at frequency /, the 
index n denotes the frame number. The quantized LSFslsf„(k)of the current frame are 
based on a 4 th order MA predcition and is given by isf „ = isf „ + a? , where lsf„ is the 
5 predicted LSFs of the current frame (a function of ^",a" 2 ,a"„a^}), and a" is the 
quantized prediction error at the current frame. The prediction error is given by 
G Af = lsf„ -lif „■ In one embodiment, the prediction error from the 4 th order MA prediction 
111 is quantized with three ten (10) dimensional codebooks of sizes 7 bits, 7 bits, and 6 bits, 
f: respectively. The remaining bit is used to specify either of two sets of predictor 
JO coefficients, where the weaker predictor improves or reduces error propagation during 
W channel errors. The prediction matrix is fully populated. In other words, prediction in 

l -i i 

p{ both time and frequency is applied. Closed loop delayed decision is used to select the 

predictor and the final entry from each stage based on a subset of candidates. The 

number of candidates from each stage is ten (10), resulting in the future consideration of 

15 10, 10 and 1 candidates after the 1 st , 2 nd , and 3 rd codebook, respectively. 

After reconstruction of the quantized LSF vector as described above, the ordering 

property is checked. If two or more pairs are flipped, the LSF vector is declared erased, 

and instead, the LSF vector is reconstructed using the frame erasure concealment of the 

decoder. This facilitates the addition of an error check at the decoder, based on the LSF 

20 ordering while maintaining bit-exactness between encoder and decoder during error free 

conditions. This encoder-decoder synchronized LSF erasure concealment improves 
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performance during error conditions while not degrading performance in error free 
conditions. Moreover, a minimum spacing of 50Hz between adjacent LSF coefficients is 
enforced. 

As shown in FIG. 2, the pre-processed speech 207 further passes through a 

5 perceptual weighting filter module 228. According to one embodiment of the present 

invention, the perceptual weighting filter module 228 includes a pole zero filter and an 

adaptive low pass filter. The traditional pole-zero filter is derived from the unquantized 

O A( y ) 

LPC filter given by: WAz)= /r ) , where y x = 0.9 and y 2 =0.55. The pole-zero filter is 

primarily used for the adaptive and fixed codebook searches and gain quantization, 
fo The adaptive low-pass filter of the module 228, however, is given by 

m w 2 (z) = — 5 — , where 77 is a function of the tilt of the spectrum or the first reflection 

hi \-r]z 1 

2 coefficient of the LPC analysis. The adaptive low-pass filter is primarily used for the 
open loop pitch estimation, the waveform interpolation and the pitch pre-processing. 

Referring to FIG. 2, the encoder 200 further classifies the pre-processed speech 

15 signal 207. The calssification module 230 is used to emphasize the perceptually 
important features during encoding. According to one embodiment, the three main 
frame-based classifications are detection of unvoiced noise-like speech, a six-grade signal 
characteristic classification, and a six-grade classification to control the pitch pre- 
processing. The detection of unvoiced noise-like speech is primarily used for generating 

20 a pitch pre-processing. In one embodiment, the classification module 230 classifies each 
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frame into one of six classes according to the dominating feature of that frame. The 
classes are: (1) Silence/Background Noise, (2) Noise-Like Unvoiced Speech, (3) 
Unvoiced, (4) Onset, (5) Non- Stationary Voiced and (6) Stationary Voiced. In some 
embodiments, the classification module 230 does not initially distinguish between non- 
5 stationary and stationary voiced of classes 5 and 6, and instead, this distinction is 
performed during the pitch pre-processing, where additional information is available to 
the encoder 200. As shown, the input parameters to the classification module 230 are the 

O pre-processed speech signal 207, a pitch lag 231, a correlation 233 of the second half of 

m each frame and the VAD information 225. 

lIS Turning to FIG. 2, it is shown that the pitch lag 23 1 is estimated by an open loop 

pitch estimation module 232. For each 20ms frame, the open loop pitch lag has to be 

10 estimated for the first half and the second half of the frame. These estimations may be 

%l used for searching an adaptive code-book or for an interpolated pitch track for the pitch 
pre-processing. The open loop pitch estimation is based on the weighted speech given by 

15 s M .(z) = S(z)-w,(z)w 2 (z) , where s(z) is the pre-processed speech signal 207. Two sets of 
open loop pitch lags and pitch correlation coefficients are estimated per frame. The first 
set is centered at the second half of the frame and the second set is centered at the first 
half frame of the subsequent frame, i.e. the look-ahead frame. The set centered at the 
look-ahead portion is recycled for the subsequent frame and used as a set centered at the 

20 first half of the frame. Accordingly, for each frame, there are three sets of pitch lags and 
pitch correlation coefficients available to the encoder 200 at the computational expense of 
only two sets, i.e., the sets centered at the second half of the frame and at the look-ahead. 
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Each of these two sets is calculated according to the following normalized correlation 

i 

function: /?(&) = ^ , where z = 80 is the window size, and £ = ^s w («) 2 is the 

energy of the segment. The maximum of the normalized correlation R(k) in each of three 
regions [17,33], [34,67], and [68,127] are determined, which determination results in 
5 three candidates for the pitch lag. An initial best candidate from the three candidates is 
selected based on the normalized correlation, classification information and the history of 
p the pitch lag. 

* ; ! Turning back to the speech pre-processor block 210, as discussed above, the noise 

suppression module 206 receives various voicing parameters from the speech processor 

W) block 250 in order to improve the calcualtion of the channel gain. The voicing 
parametres may be derived from various modules within the speech processor block 250, 

O such as a the classification module 230, the pitch estimation module 232, etc. The noise 

^ suppression module 206 uses the voicing parametrs to adjust the channel gains {y ch (i)} ■ 

As explained above, the goal of noise suppression, for a given channel, is to adjust 

15 the gain y ch such that it is higher or closer to 1.0 to preserve the speech quality for strong 
voiced areas and, on the other hand, lowring the gain y ch to be closer to zero for 
suppressing the noise in noisy areas of speech. Theoratically, for a pure voice signal, the 
gain y ch should be set to "1.0", so the signal remains unmodifed, on the other hand, for a 
pure noise signal, the gain y ch should be set to "0", so the noise signal is suppressed. In 

20 between these two theoratical extrems, there lies a spectrum of possible gains y ch , where 
for voice signals, it is desirable to have a gain y ch closer to "1.0" to preserve the speech 
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quality as much as possible. Now, since the speech processor block 250 contributes to 
cleaning or suppressing some of the noise in the voiced areas, the conventional noise 
suppression process may be relaxed (as discussed below.) For example, referring to FIG. 
3, speech sections 302, 304 and 306 that are located between the harmonics in the voiced 
5 area have a very low signal-to-noise ratio and as a result the speech sections 302, 304 and 
306 are noisy sections of the voiced area. But, it should be noted that the speech 
processor block 250 contributes to cleaning the noisy speech areas 302, 304 and 306 by 

IS applying pitch enhancement. Accordingly, modification of the speech signal by reducing 

: J1 the gain y ch in such areas may be avoided. 

4p The present invention overcomes the drawbacks of the conventional approaches 

f and improves the gain computation by using other dynamic or voicing parameters, in 
} addition to the SNR parameter used in conventional approaches to noise suppression. In 

t""1 

ti one embodiment of the present invention, the voicing paramters are fed back from the 
speech precossor block 250 into the noise suppression module 206. These voicing 

15 parametrs belong to previously processed speech frame(s). The advantage of such 
embodiment is achieving a less complex system, since such embodiment reuses the 
information gathered by the speech processor block 250. In other embodiments, however, 
the voicing parametrs may be calculated within the noise suppression module 206. In 
such embodiments, the voicing parameters may belong to the particular speech frame 

20 being processed as well as those of the preceding speech frames. 

Regardless of whether the voicing parameters are fed back to the noise suppression 
module 206 or are calculated by the noise suppression module 206, in one embodiment, 
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the channel gain is first calculated in the db domain based on the following equation: 
7db(i) = MgO) (v'VO - °th) + 7m where the gain slope // g (i) is defined as: 



Mg<J) = { 



0.39, if voicing parameters indicate unvoiced speech 

0.39+x , where 0 < x < 0.61, if voicing parameters 
indicate voiced speech 



Yet, in other embodiments, the voicing parameters may be used to modify any of 
the other parameters in the y db (i) equation, such as y n or o lh . Nevertheless, the voicing 
10 parameters are used to adjust the gain for each channel through the calculation of the 
3 value of "x" by the noise suppression module 206. For example, in one embodiment, the 
F noise suppression module 206 may use the classification parameters from the 
C classification module 230 to calculate the adjustment value "x". As explained above, in 
L one embodiment, the classification module 230 classifies each speech frame into one of 
jfi5 the six classes in accordance to the dominating features of each frame. With reference to 
C FIG. 4, if the frames is classified to be in the unvoiced area 410, ju g Q) will be °- 39 - 
However, if the frames is classified as being in the voiced area 420, ju^i) will 0.39 + x, 
and "x" may be adjusted based on the strength of the voice signal. For example, if the 
voice signal is classified as stationary voiced, the value of "x" will be higher, but for non- 
20 stationary voiced classification, the value of "x" will be less. 

In addition to the classification parameter, one embodiment may also consider the 
pitch correlation R(k) . For example, in the voiced area 420, if the pitch correlation value 
is higher than average, the value of "x" will be increased, and as a result the value of fi g (i) 
is increased and the speech signal G(k) is less modified. Furthermore, an additional 
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factor to consider may be the value of // g (i-l), since the value of /u g (i) should not be 
dramatically different than the value of its preceding n z . 

The present invention may be embodied in other specific forms without departing 
from its spirit or essential characteristics. The described embodiments are to be 
considered in all respects only as illustrative and not restrictive. For example, the voicing 
parameters that are calculated in the speech processing block 250 may be used or 
considered in a variety of ways and methods by the noise suppression module 206 and the 
present invention is not limited to using the voicing parameters to adjust the value of 
some parameters, such fi g , y n or a th . The scope of the invention is, therefore, indicated by 
the appended claims rather than the foregoing description. All changes which come 
within the meaning and range of equivalency of the claims are to be embraced within 
their scope. 
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CLAIMS 

1 . A method of suppressing noise in a signal, said method comprising the 
steps of: 

estimating a signal to noise ratio for said signal; 
5 classifying said signal to a classification; 

calculating a gain for said signal using said signal to noise ratio and said 
classification; and 
X modifying said signal using said gain. 

I? j 2. The method of claim 1 further comprising a step of estimating a pitch 

Ip correlation for said signal, wherein said calculating step further uses said pitch 
s correlation. 

03 3 . The method of claim 1 , wherein said signal is one channel of a plurality of 

□ channels of a speech signal. 

4. The method of claim 1, wherein said calculating step calculates said gain 
1 5 based on y db = fi g ( a" q - a th ) + y n , and wherein fi g is adjusted according to said 

classification. 

5 . The method of claim 2, wherein said calculating step calculates said gain 
based on y db = fi g (a" q -o th )+ y„, and wherein n % is adjusted according to said 
classification and said pitch correlation. 

20 6. The method of claim 1 , wherein said signal is in a time domain, and said 

method further comprises a step of converting said signal from said time domain to a 
frequency time prior to said estimating step. 
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7. The method of claim 1, wherein said signal is in a frequency domain, and 
said method further comprising a step of converting said signal from said frequncy 
domain to a time domain after said modifying step. 

8. A method of suppressing noise in a signal having a first portion and a 
second portion, said method comprising the steps of: 

computing a voicing parameter using said first portion; 
estimating a signal to noise ratio for said second portion; 

calculating a gain for said second portion using said signal to noise ratio and said 
voicing parameter; and 

modifying said signal using said gain. 

9. The method of claim 8, wherein said first portion of said signal is ahead of 
said second portion in a time domain. 

10. The method of claim 8 5 wherein said voicing parameter is computed by a 
speech coder. 

1 1 . The method of claim 8, wherein said voicing parameter is a speech 
classification of said first portion. 

12. The method of claim 8, wherein said voicing parameter is a pitch 
correlation of said first portion. 

13. The method of claim 8, wherein said calculating step calculates said gain 
based on y db = ju g ( o" q - o th ) + y m and wherein ju g is adjusted according to said voicing 
parameter. 

14. The method of claim 8, wherein said calculating step calculates said gain 
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based on y db = ju g ( c" q - o th ) + y m and wherein y n is adjusted according to said voicing 
parameter. 

15. The method of claim 8, wherein said signal is in a time domain, and said 
method further comprises a step of converting said signal from said time domain to a 
frequency time prior to said estimating step. 

16. The method of claim 8, wherein said signal is in a frequency domain, and 
said method further comprising a step of converting said signal from said frequncy 
domain to a time domain after said modifying step. 

17. A noise suppression system comprising: 
a signal to noise ratio estimator; 

a signal classifier; 

a signal gain calculator; and 

a signal modifier; 

wherein said estimator estimates a signal to noise ratio of said signal, said signal is 
given a classification using said signal classifier, said signal gain is calculated based on 
said signal to noise ratio and said classification using said calculator, and wherein said 
signal modifier modifies said signal by applying said gain. 

18. The system of claim 17 further comprising a signal pitch estimator for 
estimating a pitch correlation of said signal for use by said gain calculator. 

19. The system of claim 17 further comprising a frequency-to-time converter to 
convert said signal from a frequency domain to a time domain. 

20. A system capable of suppressing noise in a signal having a first portion and 
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a second portion, said system comprising: 

a signal processing module for computing a voicing parameter of said first portion; 
a signal to noise ratio estimator; 
a signal gain calculator; and 
a signal modifier; 

wherein said estimator estimates a signal to noise ratio of said second portion, said 
second portion gain is calculated based on said signal to noise ratio and said voicing 
parameter using said calculator, and wherein said signal modifier modifies said second 
portion by applying said gain. 

2 1 . The system of claim 20, wherein said first portion of said signal is ahead of 
said second portion in a time domain. 

22. The system of claim 20, wherein said signal processing module is a speech 

coder. 

23. The system of claim 20, wherein said voicing parameter is a speech 
classification of said first portion. 

24. The system of claim 20, wherein said voicing parameter is a pitch 
correlation of said first portion. 

25. The system of claim 20, wherein said signal gain calculator determines said 
gain based on y db = p g (c" q - <J th ) + y m and wherein ju g is adjusted according to said 
voicing parameter. 

26. The system of claim 20 further comprising a frequency-to-time converter to 
convert said second portion of said signal from a frequency domain to a time domain. 
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ABSTRACT 

An input signal enters a noise suppression system in a time domain and is 
converted to a frequency domain. The noise suppression system then estimates a signal to 
noise ratio of the frequency domain signal. Next, a signal gain is calculated based on the 
estimated signal to noise ratio and a voicing parameter. The voicing parameter may be 
determined based on the frequency domain signal or may be determined based on a signal 
ahead of the frequency domain signal with respect to time. In that event, the voicing 
parameter is fed back to the noise suppression system, for example, by a speech coder, to 
calculate the signal gain. After calculating the gain, the noise suppression system 
modifies the signal using the calculated gain to enhance the signal quality. The modified 
signal may further be converted from the frequency domain back to the time domain for 
speech coding. 



-30- 



00CXT0584C 



105 



104 



High-Pass Filter 



103 



102 



Silence 
Enhancement 



Input Speech 



101 



Noise Suppression 



Frequency 

Domain 
Conversion 



Spectral 
Deviation . 
Estimata 



I 



-MX 



r 1 




Channel 
Energy 
Estimatcr 



Channel 
Gain 
Calculation 



^0 



Time 

Domain 
Conversion 



j± 



Channel ^Ajr\ x 

SNR 
Estimator 



Vdce 
Metric 
Calculation 



Sin) 



Background 

Noise 
Estimator 



\H0 



Channel 

SNR 
Modifier 



r 



Noise 
Update 
Decision 



— i updatejlag 



FIG. 1a 
(PRIOR ART) 



/* Normal update logic */ 

update Jlag « FALSE 

if ( v(m) < UPDATE_THLD ) { 

updatejlag = TRUE 

update_cnt ~ 0 

} 

/* Forced update logic */ 

else if (( E M (m) > NOISE_FLOOR_DB ) and ( A^m) < DEV_THLD )) { 
update_cnl = update_cnt + 1 
if ( update_cnl > UPDATE_CNT_THLD ) 
update Jlag = TRUE 

} 

/* "Hysteresis" logic to prevent long-term creeping of update„cnt */ 
if ( update_cnt ~ last_update_cnt ) 

hyster_cnt = hyster_cnt + 1 
else 

hyster_cnt = 0 
lastjupdate^cnt = update_cnt 
if ( hyster_cnt > HYSTER_CNTTHLD ) 

updatejcnt = 0 



FIG. 1b 
(PRIOR ART) 



/* Set or reset modify flag */ 
index_cnt = 0 

for ( i = N u to N c - 1 step 1 ) { 
if (a,(i)>INDEX_THLD) 
indexjcni = index_cnt + 1 

} 

if ( index_cnt < INDEX_CNT_THLD ) 

modify Jlag = TRUE 
else 

modify Jlag = FALSE 

/* Modify the SNR indices to get \a' q }*/ 

if ( modify Jlag == TRUE ) 
for ( i = 0 to N c - 1 step 1 ) 

if (( v(m) < METRIC.THLD ) or (a,(i) < SETBACK.THLD )) 
1 

else 

else 

/* Limit jcr*} to SNR threshold a !h */ 
for ( i = 0 to 7V C - 1 step 1 ) 

else 

^(i)=a-,(i) 

FIG. 1c 
(PRIOR ART) 
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United States Patent Application 
COMBINED DECLARATION AND POWER OF ATTORNEY 

As a below named inventor I hereby declare that: my residence, post office address and citizenship are as stated below next to my name; that 

I verily believe I am the original, first and sole inventor (if only one name is listed below) or a joint inventor (if plural inventors are named below) of 
the subject matter which is claimed and for which a patent is sought on the invention entitled: FREQUNCY DOMAIN NOISE SUPPRESSOR 

The specification of which 

a. XX is attached hereto ^ ^ w . At _ . 

b was filed on as application serial no. and was amended on (if applicable) (in the case of a PC I- 

filed application) described and claimed in international no. filed and as amended on (if any), which 

1 have reviewed and for which I solicit a United States patent. 

I hereby state that I have reviewed and understand the contents of the above-identified specification, including the claims, as amended by any 
amendment referred to above. 

I acknowledge the duty to disclose information which is material to the examination of this application in accordance with Title 37, Code of Federal 
Regulations, Section 1.56 (see the last page attached hereto). 

I hereof claim foreign priority benefits under Title 35, United States Code, Sections 1 19/365 of any foreign application(s) for patent or inventor's 
certific&le listed below and have also identified below any foreign application for patent or inventor's certificate having a filing date before that of the 
appliclflon on the basis of which priority is claimed: 

a. XX r no such applications has been filed. 

b. [ |uch applications have been filed as follows: 



^ COUNTRY 


APPLICATION NUMBER 


DATE OF FILING 
(day, month, year) 


DATE OF ISSUE 
(day, month, year) 


















ALLlpOREIGN APPLICATIONS, IF 


- ANY. FILED BEFORE THE PRIORITY APPLICATION(S) 


COUNTRY 


APPLICATION NUMBER 


DATE OF FILING 
(day, month, year) 


DATE OF ISSUE 
(day, month, year) 



















I hereby claim the benefit under Title 35, United States Code, Sections 120/365 of any United States and PCT international application® (listed 
below and insofar as the subject matter of each of the claims of this application is not disclosed in the prior United States application m the 
manner provided by the first paragraph of Title 35, United States Code, Section 1 1 2, I acknowledge the duty to disclose material information as 
defined in Title 37, Code of Federal Regulations, Section 1 .56(a) which occurred between the filing date of the prior application and the national or 
PCT international filing date of this application. 



U.S. APPLICATION NUMBER 


DATE OF FILING (day, month, year) 


STATUS (patented, pending, abandoned) 
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I hereby appoint the following attorney(s) and/or patent agent(s) to prosecute this application and to transact ail business in the Patent and 
Trademark Office connected herewith: 



MICHAEL FARJAMI, Reg. No. 38,135 
FARSHAD FARJAMI, Reg. No. 41,014 
DANIEL N. YANNUZZI, Reg. No. 36,727 
SEMION TALPALATSKY, Reg. No. 35,380 



JAMES K. DAWSON, Reg. No. 41,701 
KELLY H. HALE, Reg. No. 36,542 
ROBERT P. HART, Reg. No. 35,184 
KEITH KIND, Reg. No. 42,735 



I hereby authorize them to act and rely on instructions from and communicate directly with the person/assignee/firm/organization who/which first 
sends/sent this case to them and by whom/which I hereby declare that I have consented after full disclosure to be represented unless/until I 
instruct them to the contrary. 

Please direct all correspondence in this case to FARJAMI & FARJAMI LLP at the address indicated below: 

FARJAMI & FARJAMI LLP 
16148 Sand Canyon 
Irvine, California 92618 
Telephone: (949) 784-4600 

I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are believed 
to be true; and further that these statements were made with the knowledge that willful false statements and the like so made are punishable by 
fine oc imprisonment, or both, under Section 1 001 of Title 18 of the United States Code and that such willful false statements may jeopardize the 
validity if the application or any patent issued thereon. 
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37 CF.R. Section L56 - Duty to disclose information material to patentability. 

A patent by its very nature is affected with a public interest. The public interest is best served, and the most 
effective patent examination occurs when, at the time an application is being examined, the Office is aware of 
and evaluates the teachings of all information material to patentability. Each individual associated with the 
filing and prosecution of a patent application has a duty of candor and good faith in dealing with the Office, 
which includes a duty to disclose to the Office all information known to that individual to be material to 
patentability as defined in this section. The duty to disclose information exists with respect to each pending 
claim until the claim is cancelled or withdrawn from consideration, or the application becomes abandoned. 
Information material to the patentability of a claim that is cancelled or withdrawn from consideration need not 
be submitted if the information is not material to the patentability of any claim remaining under consideration 
in the application. There is no duty to submit information which is not material to the patentability of any 
existing claim. The duty to disclose all information known to be material to patentability is deemed to be 
satisfied if all information known to be material to patentability of any claim issued in a patent was cited by the 
Office or submitted to the Office in the manner prescribed by Sections 1.97(b)-(d) and 1.98. However, no 
patent will be granted on an application in connection with which fraud on the Office was practiced or 
attempted or the duty of disclosure was violated through bad faith or intentional misconduct. The Office 
encourages applicants to carefully examine: 

Prior art cited in search reports of a foreign patent office in a counterpart application, and 

asa The closest information over which individuals associated with the filing or prosecution of a patent application 
H believe any pending claim patentably defines, to make sure that any material information contained therein is 
^ disclosed to the Office. 

I-:! Under this section, information is material to patentability when it is not cumulative to information already of 
f-« record or being made of record in the application, and 

■sj It establishes, by itself or in combination with other information, a prima facie case of unpatentability of a 
ff[ claim; or 

It refutes, or is inconsistent with, a position the applicant takes in: 
ft Opposing an argument of unpatentability relied on by the Office, or 
J;f Asserting an argument of patentability. 

D A prima facie case of unpatentability is established when the information compels a conclusion that a claim is 
unpatentable under the preponderance of evidence, burden-of-proof standard, giving each term in the claim its 
broadest reasonable construction consistent with the specification, and before any consideration is given to 
evidence which may be submitted in an attempt to establish a contrary conclusion of patentability. 

Individuals associated with the filing or prosecution of a patent application within the meaning of this section 
are: 

Each inventor named in the application; 

Each attorney or agent who prepares or prosecutes the application; and 

Every other person who is substantively involved in the preparation or prosecution of the application and who 
is associated with the inventor, with the assignee or with anyone to whom there is an obligation to assign the 
application. 

Individuals other than the attorney, agent or inventor may comply with this section by disclosing information to 
the attorney, agent, or inventor. 
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